SOFIA FORCAST FAR-IR PHOTOMETRY OF COMET ISON AND CONSTRAINTS ON THE COMA 
GRAIN SIZE DISTRIBUTION. D. H. Wooden 1 , J. M. De Buizer 2 , M. S. Kelley 3 , C. E. Woodward 4 , D. E. 
Harker 5 , W. T. Reach 2 , M. L. Sitko 6 , R. W. Russell 7 , R. D. Gehrz 4 , Imke de Pater 8 , Ludmilla Kolokolova 3 

! NASA Ames Research Center, Moffett Field, CA 94035, 2 SOFIA, USRA, 3 University of Maryland, College Park, 
MD, 4 MN Institute of Astrophysics, University of Minnesota, Minneapolis, MN, 5 CASS/UCSD, San Diego, CA, 
6 Space Science Institute, Cinncinatti, OH, 7 The Aerospace Corporation, El Segundo, CA, 8 University of California 
Berlekey, Berkeley, CA 


Introduction: Comet C/2012 SI (ISON) was unique 
in that it was a dynamically new comet derived from 
the nearly isotropic Oort cloud reservoir of comets 
with a sun-grazing orbit. Infrared (IR) observations 
were executed on NASA's Stratospheric Observatory 
For Infrared Astronomy (SOFIA) by the FORCAST [1] 
instrument on 2013 October 25 UT (r h =1.18 AU, 

A=1 .5AU). Photometry was obtained in FORCAST 
filters centered at 11.1, 19.7, and 31.5 pm. The obser- 
vations compliment a large world-wide effort to ob- 
serve and characterize comet ISON. 


Comet dust properties, including the dust mineralogy, 
porosity, and grain size distributon parameters are de- 
rived by fitting thermal dust models to IR observations 
of comets. Presently, the dust properties are well con- 
strained for relatively few long-period Nearly Isotropic 
Comets (NICs) dynamically from the Oort cloud com- 
pared to short-period Ecliptic Comets (ECs) from the 
scattered disk. Previously studied NICs have dust 
properties that span a range of small and/or highly po- 
rous grains (e.g., C/1995 01 (Hale-Bopp) [3,4,5] and 
C/2001 Q4 (NEAT) [6] ) to large and/or compact grains 
(e.g., C/2007 N4 (Lulin) [7] and C/2006 PI (McNaught) 
[8] ). For the sample of NICs, silicate crystalline mass 
fractions range from 7*^-0. 1-0.6 [8,9] , which are 
benchmarks for radial transport in our protoplanetary 
disk [10] . Comet Hale-Bopp near perihelion revealed 
very strong silicate features, with a silicate strength of 
>3 and had crystalline peaks near 16, 19, 23.5, 27.5, 

33 pm, [11] which yield a high crystalline fraction of 
/ayst~0.6. [3] Furthermore, Hale-Bopp’ s forsterite fea- 
tures’ central wavelength positions, relative intensities, 
and feature asymmetries can constrain the shapes of 
the forsterite crystals; the crystal shapes are consistent 
with high temperature disk condensates and evapo- 
rates. [12] The grain size distribution parameters as well 



Figure 1. SOFIA FORCAST imaging photometry of 
comet ISON in 3 filters: 11.1, 19.7, and 31.5 pm. 


as the mineralogy give clues to comet origins. In ra- 


dial transport models for dust grains in the protoplane- 
tary disk, larger grains (~20 pm) are harder to transport 
to larger disk distances (>~20 AU) than smaller grains 
(~1 |am). [13] 

Dust Modeling Methods: Through modeling of ob- 
servations of 1 0 pm silicate emission feature and its 
strength, [2] the mineralogy can be constrained. For ob- 
servations with sufficient signal-to-noise ratio (SNR), 
models can constrain the silicate-to-carbon ratio and 
the olivine-to-pyroxene ratio. If a 1 1.1-1 1.2 pm peak 
from forsterite is detected, then the models further con- 
strain the silicate crystalline-to-amorphous ratio, or 
crystalline fraction / cryst = m cryst /(m amor + m cryst ). 

Through modeling the mid- and far-IR flux densities, 
we constrain the properties of the grain size distribu- 
tion (GSD) and grain porosity (P versus dust radius a, 
parameterized by D such that P = (a/0.1 pm) A (D-3), 

D= 3 for solid and D= 2.5 for highly porous). [7] In a 
Hanner (modified power-law) grain size distribution, [7] 
the small radii grains at the peak of the grain size dis- 
tribution (tf p ) dominate the surface area and the flux 
density. Smaller grains produce higher contrast silicate 
features. Grains of greater porosity grains also produce 
higher contrast silicate features. Similarly, larger 
grains are cooler, contribute to the ‘continuum’ under 
the silicate feature and contribute to the far-IR flux 
density. A steeper grain size distribution has more 
smaller grains relative to larger grains, so a steeper 
grain size distribution can produce a stronger silicate 
feature because the smallest grains dominate the dust 
emitting surface area. More smaller grains also can 
produce a stronger ~3. 5-8 pm ‘featureless’ continuum, 
which is typically attributed to submicron amorphous 
carbon grains. [2,8,9,10,14] The amorphous carbon feature- 
less emission extends through the 10 pm region, so a 
higher amorphous carbon-to-silicate ratio also can 
weaken the silicate feature. [6] As described by these 
examples, the mineralogy and grain size distribution 
parameters (porosity, a p , and slope N) are coupled and 
must be sought simultaneously when fitting thermal 
models to comet IR observations. Even with these 
modeling complexities, given measurements of the 
silicate feature strength, the color temperature of the 
10 pm continuum, and three photometry measurements 
near 11.1, 19.7, and 31.5 pm, this is sufficient informa- 





tion to constrain the silicate-to-carbon ratio and grain 
size distribution parameters (porosity, a v , and slope TV). 

Dust Models of Comet ISON (r h ~1.2 AU): 

For modeling comet ISON at this epoch, we consider 
three mineral compositions: amorphous carbon, and 
amorphous silicates with 50% olivine and 50% pyrox- 
ene. We consider a grain size distribution (GSD) 
spanning grain radii 0.1 <a< 1000 pm, slopes 3.2 < TV 
< 4.5 in steps of 0.1, and four porosities characterized 
by porosity parameter [7] D= 2.500, 2.727, 2.857, and 
3.000. When the thermal models are fitted to the 
FORCAST photometry points and employ the con- 
straint on the silicate feature strength (<1.1) given from 
ground-based N-band spectra by T. Ootsubo and col- 
laborators, [15] then lower reduced goodness-of-fit (x v 2 ) 
values are achieved for smaller a v and steeper slopes 
N, as shown in the four-panel Fig. 2. If the constraint 
of color temperature of the N-band continuum is added 
to the thermal model fits, then low values are 
achieved for larger a v (0.7pm < a p < 1.0pm), porosity 
parameter P=2. 7, and an even steeper GSD slope 
(TV=4.5), as shown in Fig. 3. The N-band color tem- 
perature is an important constraint on the GSD. 



Figure 2: Reduced chisq (y v 2 ) for thermal models fit- 
ted to three FORCAST photometry points and with the 
constraint that the silicate feature strength <1.1 for 
comet ISON at r b ~1.2 A U. 

Results: Thermal model fits to SOFIA +FORCAST 
broadband photometry of comet ISON at three wave- 
lengths, 11.1, 19.7 and 31.5 pm, offer constraints on 
the coma grain size distribution when combined with a 
limit on the silicate feature strength (<1.1) and color 
temperature (260-265 K) of the continuum on the short 
and long wavelength sides the silicate feature, as 
measured within a week of FORCAST photometry by 
Subaru COMICS ground-based N-band spectros- 
copy. [15] For this epoch (r h ~1.2 AU), comet ISON is 
better fitted by thermal modesl with dust properties 
including moderately porous grains (D=2.7), a peak 
grain radius of 0.7pm < a v <1 pm, and a size distribu- 


tion slope of TV=4.5. The peak grain radius of comet 
ISON at this epoch is more similar to two Oort clouds 
comets with large/and or compact porous grains 
(C/2007 N4 (Lulin) m and C/2006 PI (McNaught) [8] ), 
but its size distribution slope is considerably steeper 
than all other NICs. For such a steep grain size distri- 
bution to generate such a weak silicate feature requires 
little silicates compared to carbon. Specifically, amor- 
phous carbon : silicate - 0.9 : 0.1. At this epoch, 
comet ISON’s amorphous carbon-to-silicate ratio is 
high compared to other NICs from the Oort cloud and 
is on the high end of the range of all comets including 
Ecliptic Comets from the scattered disk [8] . Comet 
ISON’s thermal emission is domainted by -0.7-1 pm 
grains that are carbon rich. 
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Figure 3: Better-fit thermal models for comet ISON at 
r h ~1.2 AU, for a p =0. 7pm (left) and a p =1.0gm (right), 
porosity parameter D=2 . 7 and slope N=4.5. Models 
are fitted to three FORCAST photometry points (black 
points with error bars, and constrained by the silicate 
feature strength <1.1 and the N-band continuum color 
temperature ofT bb =260-265 A^ l5] (blue dashed line of 
a graybody at this temperature). Mineralogy is amor- 
phous carbon-to-amorphous silicate ratio = 9, with 
amorphous carbon (orange), amorphous pyroxene, 
amorphous olivine (green), and the total SED (red). 
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